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ABSTRACT 


The  problem  considered  is  that  of  determining  the  radius  of 

oo 

convergence  of  the  series  S  u(r)x  ,  when  the  coefficients  u(r)  are 

0 

connected  by  a  recurrence  relation  containing  (p  +  1)  terms.  Apart 
from  the  trivial  case  P  =  1,  only  the  case  p  =  2  has  hitherto  been 
treated,  using  the  method  of  continued  fractions,  which  cannot  be 
generalized  to  p  >  2. 

In  this  paper,  a  matrix  technique  is  given,  which  is  applicable 
to  any  value  of  £  ,  whereby  infinite  products  of  (p  X  p)  matrices  are 
used  to  formulate  characteristic  equations  and  expressions  for  the  u(r). 
Particular  attention  is  given  to  the  situation  in  which  the  problem 
involves  (p  -  1)  variable  parameters,  whose  values  have  to  be  deter¬ 
mined  in  order  to  secure  the  greatest  possible  radius  of  convergence. 


CONVERGENCE  OF  RECURRING  SERIES  WITH  GENERAL-ORDER 


RECURRENCE  RELATIONS 
F.  M.  Arscott 


1.  Introduction 

In  this  paper  we  investigate  convergence  properties  of  a  series 

,  (1.1) 

0 

when  there  is  a  recurrence  relation  between  (p  +  1)  successive 
coefficients  u(r)  of  the  series,  of  the  form 

u(r+ 1)  =  aj(r)u(r)  +  0-2 ®pWu(r- p+ 1)  ,  (r>p-l)  (1.2) 

in  which  the  a^(r)  have  known  asymptotic  behavior  as  r  -►  oo  ,  The  main 

feature  of  interest  is  the  phenomenon  which  I  shall  call  "augmented 

convergence"  ;  in  general,  the  series  (1. 1)  will  have  a  certain  radius  of 

convergence  Pj  ;  if  a  certain  relation  holds  between  the  coefficients 

a^(r)  and  the  "starting  values"  u(0),  u(l), .  .  .  ,u(p  -  1),  then  the  radius 

of  convergence  may  be  Increased  to  a  higher  value  p  ;  if  a  further 

relation  is  also  satisfied  then  the  radius  of  convergence  may  be  Increased 

to  a  still  higher  value  p^,  and  so  on.  Ultimately,  if  (p  -  1)  relations 

are  satisfied,  the  radius  of  convergence  reaches  its  greatest  possible 

value  p  . 

*P 

Sponsored  by  the  Mathematics  Research  Center,  United  States  Army, 
Madison,  Wisconsin  under  Contract  No.  DA-l  1-022-ORD-2059. 


-2- 


#349 


The  investigation  of  this  problem  was  stimulated  by  the  occurrence 
of  series  such  as  (i.  i)  In  the  solution  of  ordinary  differential  equations 
by  the  method  of  Frobenius.  In  the  classical  equations  of  hypergeometric 
type  (Legendre's  and  Bessel's  equations,  for  Instance),  p  =  1  so  the 
problem  is  trivial;  D'Alembert's  ratio  test  gives  the  radius  of  convergence 
unambiguous  as  |lim  o^Cr)  |  ^  ,  and  the  phenomenon  of  augmented  con¬ 
vergence  does  not  occur.  But  in  equations  of  higher  type  than  hypergeo¬ 
metric  (Mathieu' s  and  Lame' s  equations,  the  spheroidal  and  ellipsoidal  wave 
equations,  etc.)  the  recurrence  relations  always  Involve  at  least  three 
terms,  often  four  or  five.  Then,  the  securing  of  augmented  convergence 
is  equivalent  to  choosing  the  parameters  in  the  differential  equation  so 
that  a  solution  may  converge  in  a  region  containing  singularities  of  the 
equation. 

Two  forms  in  which  the  problem  commonly  arises  are 

(i)  the  a^(r)  of  (1. 2)  are  all  fully  determined  functions  of  jr  ,  while 
the  "starting  values"  u(0), . . . ,  u(p  -  1)  are  (initially  at  least)  arbitrary; 

(ii)  the  starting  values  are  given  in  terms  of  a  single  value  u(0) 
by  (p  -  1)  supplementary  relations  of  the  form 

■u(l)  =  aj(0)u(0) 

u(2)  =  a,(l)u(l)  +  a_(l)u(0)  , 

(1.2a) 

u(p  -  1)  =  o  (p  -  2)u(p  -  2)  +  . . .  +  Of  (p  -  2)u(0)  , 

1  P“  1  ' 

and  the  functions  o^(r)  contain  (p  -  1)  parameters.  Each  of  these  special 
cases,  and  the  intermediate  situation  in  which,  say,  the  functions  a^(r) 
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contain  (p  -  k)  parameters  and  there  are  (k  -  1)  relations  of  the  form 
(1.2a),  are  Included  in  what  follows. 

Examples  of  three-term  recurrence  relations  (the  case  p  =  2) 
occurred  in  the  latter  part  of  the  19th  century.  These  were  treated  by 
Kelvin  and  Heine  (apparently  independently)  using  a  continued-fraction 
technique  which  has  been  extensively  developed  by  Perron  (9,  Vol.II  §  20) 
Meixner  and  Schafke  (7,  §  1.8)  Blanch  (5),  Bouwkamp  (6)  and  others  to 
the  point  where  it  provides  a  completely  satisfactory  treatment  of  the 
case  p  =  2,  both  from  the  theoretical  and  the  computational  points  of  view 

The  only  drawback  to  the  continued-fraction  technique  is  the 
apparent  impossibility  of  direct  generalization  to  values  of  £  greater  than 
2.  This  is  hardly  surprising,  since  the  essential  feature  of  a  continued- 
fraction  is  that  the  successive  convegents  have  the  form  >  where 

and  B^  each  satisfy  recurrence  relations  with  precisely  three  terms. 
Moreover,  a  careful  study  of  the  technique  suggests  that  the  use  of 
continued  fractions  provides  a  convenient  formulation  of  the  problem  but 
little  positive  help  towards  a  solution,  since  in  order  to  handle  such 
fractions  one  is  thrown  back  ultimately  on  to  a  (possibly  different) 
three-term  recurrence  relation. 

In  this  paper,  therefore,  the  continued-fraction  method  has  been 
abandoned;  instead  we  follow  the  lead  given  thirty  years  ago  by 
Milne-Thomson  (8,  §  5.  3)  and  by  George  Blrkhoff  and  his  co-workers  (1-4) 
and  use  a  matrix  formulation  of  the  problem.  This  provides  an  alternative 
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method  of  treating  the  three-term  recurrence  relation  and  one  which 
may  well,  in  view  of  the  ease  with  which  modem  computers  handle 
matrix  multiplication,  prove  more  suitable  for  computational  purposes. 

It  also  generalizes,  very  easily  and  naturally,  to  the  case  of  recurrence 
relations  of  any  finite  order. 

The  problem  as  a  whole  has  three  inter-related  aspect s:- 
(i)  determination  of  the  possible  radii  of  convergence  of  the 
series, 

(il)  formulation  of  the  conditions  to  be  satisfied  in  order  to  secure 
augmented  convergence  (the  characteristic  equations). 

(iii)  determination  of  the  associated  values  of  the  coefficients  u(r). 
Sometimes  only  (i) ,  or  (l)  and  (li)  are  required. 

In  the  present  paper,  we  give  an  outline  of  the  general  method  in 
§  2,  then  in  §  3  re-formulate  the  problem  in  matrix  notation  and  obtain 
several  different  forms  of  the  conditions  for  augmented  convergence.  We 
are  then  able  to  justify,  in  §  4,  a  working  rule  by  which  the  possible 
radii  of  convergence  may  be  written  down  almost  Immediately.  Finally, 
in  S§  5,6  the  evaluation  of  the  coefficients  u(r)  is  considered,  a  process 
which  yields  further  forms  of  the  conditions  for  augmented  convergence. 

Besides  my  obvious  indebtedness  to  the  authors  cited,  my  wannest 
thanks  are  also  due  to  Professor  Tomlinson  Fort  for  drawing  my  attention 
to  the  relevance  of  Birkhoff's  work  to  the  problem. 

We  shall  have  occasion  to  use  the  following  results  relating  to 
difference  equations  of  type  (1.2),  due  to  Birkhoff  (2-3),  Birkhoff  and 
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Trjitzlnsky  (4)  and  Adams  (1).  Assume  that  each  coefficient  o'^(r)  can 
be  represented  asymptotically,  for  large  j; ,  by  a  series  containing  at 
most  a  finite  number  of  positive  powers  of  ^  ,  1.  e.  of  the  form 

00 

a^(r)~S  ,  (1.3) 

j=-k 

with  real  •  Then  for  sufficiently  large  j[  there  exist  £  linearly 
independent  solutions  of  (1.2)  say,  w,(r),  w_(r),  ...,w  (r). 

The  actual  determination  of  the  w^(r)  is  in  the  form  of  series 
which  are  asymptotically  valid  for  large  j[.  In  general,  each  solution 
w^(r)  is  of  the  form 

(T 

w^(r)  =  r V  t^"^  r  F^(r)  ,  (1.4a) 

where  F  (r)  is  a  series  asymptotic  in  _r  , 

^  id)  i  (1) 

F  (r)  1  +  “1  +  +  . . . .  (l,  4b) 

1  r 

The  constants  (i>  t,  v,  f^,  f^,  ....  may  be  found  successively  by 
substitution  of  the  formal  series  (1.4)  into  the  equation  (1.2).  The 
determination  of  |x  is  normally  simple;  it  is  real,  rational  and  in  many 
cases  is  either  zero  or  an  integer  which  may  be  found  by  inspection.  Then 
the  possible  values  of  X  appear  as  roots  of  an  algebraic  equation  of 
degree  not  more  than  £.  When  (i  and  X  have  been  evaluated,  <r  can 
be  found  with  a  little  more  difficulty,  and  one  can  then  proceed  to  the 
determination  of  f^,  f^, . . . .  ;  for  discussing  the  convergence  properties 
of  2  u^x"^  ,  however,  we  do  not  need  to  know  these  and  indeed  it  is 
often  sufficient  to  have  the  values  of  |i  and  ^  . 
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This,  as  has  been  said,  is  the  solution  "in  general",  and  here 

"in  general"  means  provided  the  value  of  ^  is  not  a  multiple  root  of 

the  equation  in  which  it  appears.  If  this  is  not  the  case,  however,  and 

we  have  an  N-fold  root,  then  N  solutions  are  obtained  which  may  take 

different  forms  according  to  circumstances  -  Indeed,  there  are  so  many 

possibilities  that  even  a  summary  would  be  inordinately  long.  The  case 

of  double  roots  should  be  mentioned,  however; if  we  have  a  value  of 

for  which  a  double  root  t^  occurs  then  there  are  normally  two 

corresponding  solutions  of  the  form 

1 

(1  r  “O’  (p.  <p- 

r  ^  t  exp  (±Yr^)r  ^[l  +  ~+  —  +  ....]  (1.5) 

2  ^ 

r  ‘■ 

y  being  a  further  real  constant  (1).  Even  this  rule,  however,  has 
exceptions  and  logarithmic  terms  may  occur  (3). 
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2.  Outline  of  the  method 

In  order  to  keep  the  exposition  within  reasonable  bounds,  we  shall 
make  here  the  simplefying  assumption  that  all  the  solutions  w^(r)  of 
(1.2)  are  of  the  form  (1.4)  or  the  form  (1.5).  As  has  been  mentioned 
above,  the  cases  thus  excluded  are  highly  exceptional,  and  it  seems  that 
extension  of  the  method  to  cover  these  would  not  present  any  intrinsic 
difficulty. 

The  solutions  w^(r)  can  clearly  be  ordered  in  such  a  way  that,  as 
r  “►<»,  we  have 

|w  (r)  I  >  |w  (r)  I  > . >  |w  (r)  I  (2.  1) 

1  C  p 

and  in  what  follows  we  assume  this  has  been  done. 

Now  the  equation  (1. 2)  is  linear  and  homogeneous,  so  any  solution 
u(r)  must  be  expressible  as  a  linear  combination  of  the  w^(r),  for 
sufficiently  great  _r  -  that  is  to  say,  there  must  be  constants  aj(  i  =  1  to  p) 
such  that 

P 

u(r)  =  S  a  w.(r)  (2.2) 

i=l 

for  jr  sufficiently  great;  the  constants  a^  will  depend  on,  and  be  uniquely 
determined  by,  the  values  taken  by  the  functions  a^{T),  ... ,  arp(r)and  the 
starting-values  u(0) ,  ... ,  u{p-l)  . 

Suppose,  first,  that  the  functions  w^(r)  are  such  that  Wj(r)/w^(r)-*0 
as  r -*•  00  for  i  <  j;  brief  consideration  of  the  forms  (1.4),  (1.5)  shows 
that  this  will  be  the  case  provided  we  exclude  the  case  =  Pj  > 
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|tj  =  |tj  I  ,  ^  ,  that  is,  we  do  not  have  two  solutions  with  the 

same  value  of  (x  and  values  of  t  which  differ  in  amplitude  but  not  in 
modulus.  Then  by  writing  (2.  2)  in  the  form 

P 

u(r)  =Wj(r)[aj  +  S  a^w^(r)/Wj(r)]  (2.3) 

2 

we  obtain  immediately  the  result  that  as  r  oo^ 

u(r)~ajWj(r)  ,  (2.4) 

r  r 

provided  a^  #  0,  so  that  the  series  Su(r)x  and  SWj(r)x  are 
equiconvergent.  Application  of  the  Cauchy  convergence  test  then  shows 
at  once  that  the  radius  of  convergence  of  S  u(r)x  is  ,  where 

Pj  =  llm  |wj(r)r^^*‘  ;  (2.  5) 

r-*  00 

from  (1.4),  (1.5)  this  clearly  has  the  value  0  if  >  0,  oo  if  <  0 

and  It j  I  ^  if  |ij  =  0. 

If,  however,  a^  =  0,  then  (2.4)  no  longer  holds;  Instead  of  (2.3) 
we  may  write 

P 

u(r)  =  W2(r)[a2  +  S  w^(r)/w2(r)]  ,  (2.6) 

giving 

u(r)~a2W2(r)  ,  (2.7) 

provided  ^2*  so  that  the  radius  of  convergence  of  £  u(r)x’^  becomes 
p^sllm  |w2(r)|  l.e.  0  if  »*2  >  **2  ^ 
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=  0.  This  is  certainly  not  less  than  and  will  in  fact  be  greater 

than  if  1*2  H-2  *  *^1  ^  *1  "  ®2 

a^  ^  0,  the  radius  of  convergence  becomes  p^  2  lim  Iw^W  I 

j* 

The  effect  of  this  is  that  the  radius  of  convergence  of  Z  u(r)x 
will  in  general  have  a  certain  value;  if  we  can  choose  the  starting-values 
u(0),  ...  ,Up  j  or  parameters  in  the  functions  a^(r)  or  both,  in  such  a  way 
as  to  make  the  constant  a ^  vanish,  then  the  radius  of  convergence  may 
be  increased.  If  we  can  also  make  a^  vanish  then  the  radius  of  con¬ 
vergence  may  be  increased  still  further,  and  so  on. 

Now  consider  the  case  excluded  above,  namely  when  two  or  more 
of  the  w^(r)  have  the  same  p.  and  values  of  t  which  are  equal  in 
modulus  but  differ  in  amplitude.  For  definiteness,  suppose  we  have 

ie,  i0- 

"*^2  '  ^1  ,  ^2  =  Te  ^  ,  and  w^W/w^W -*•  0  for  i  >  3. 

Then  as  r  -►  oo^ 


u(r) 


,*  ie  r 
VUje  * 


'1  ‘V  'z, 

r  +  a^e  r  ] 


(2.8) 


then  somewhat  tedious  consideration  of  inequalities  show  that  provided 

a^  and  a^  are  not  both  zero, 

ie.r  (T,  ie  r  <r  l/r 

[a^e  r  +  a^e  ^  r  -  1  ,  (2.9) 

giving 

1  /  *^l 

lim|u(r)r  =  limr  T  ,  (2. 10) 

r-*oo 

so  that  the  radius  of  convergence  of  Z  u(r)x’^  is  now 

Pj  2  lim  r  Itjl  ,  and  this  will  be  the  case  provided  a^  a^  do  not 
r  00 
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both  vanish.  If,  however,  =  a2  =  0,  then  the  radius  of  convergence 
becomes  =  lim  |w2{r)| 

The  same  argument  is  applied  without  difficulty  to  similar  cases. 
If,  for  instance,  we  have  w^(r)/Wj(r)  -►  0  for  i  =  2,  3, . . . ,  p  , 

w.W/w^W  -  0  for  i  =  4,  5,  ...,p,  but  1^2  =  1^3  >  1^2^  ”  ^^3^  » 
t^  ^  t^  then  the  radius  of  convergence  is  unless  a^  =  0, 


if  a^  =  0  and  either  ^2  ^  ^  ^  P4 


if  a,  =  a  =  a  =  0, 


etc. 
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3.  Formal  solution  of  the  problem 

As  remarked  above,  calculation  of  the  solutions  w^(r)  of  (1.2),  in 
full,  is  a  task  of  prohibitive  magnitude,  and  the  technique  to  be  employed 
makes  use  of  a  device  whereby  we  only  need  to  determine  the  dominant 
parts  of  the  w^(r),  a  comparatively  simple  matter.  We  write  v^(r)  for 
the  part  of  w^(r)  which  is  left  after  omitting  the  terms  which  are  o(l), 
and  y^(r)  for  the  part  of  v.(r)  left  after  omitting  the  r  term.  That 
is  to  say,  if  w^(r)  is  given  by  (1.4)  then 

Vi(r)  =  r  t^"^  r  ,  y.(r)  =  r  t.*^  ,  (3.  la) 

and  if  w^{r)  is  given  by  (1. 5), 

,  1 
v^(r)  =  r  t^*^  exp(Yr  2)r  ,  y^(r)  =  r  tj"^  exp(Yr  )  ,  (3.1b) 

so  that  as  r  -*■  <»  , 

“■i 

w^(r)  ~  v^(r)  ~  y^(r)  r  .  (3.  Ic) 


The  problem  becomes  amenable  to  treatment  only  wh^n  we  turn 
(1.2)  into  a  matrix  equation  of  the  first  order.  We  Introduce  column 
vectors  U(r),  A,  and  matrices  K(r),  W(r),  V(r),  Y(r)  as  follows:- 


U(r)  =  {u(r  +  1),  u(r),  u(r  -  1), . . . ,  u(r  -  p  +  2)}  , 


A  =  (Qj,  ®3»  ••  ’’^p}  * 

aj(r)  a^ix)  or^Cr)  .  ap(r) 

1  0  0  .  0 

0  1  0  .  0 


0  0  0  ...  1  0 


(3.  2) 
(3.  3) 


(3.4) 
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(Observe  that  det  K(r)  =  (r)  in  (3.4)) 

P 


w^(r+l) 

W2(r+1) 

..  w  (r+1) 

P 

Wj(r) 

W2(r) 

. .  Wp(r) 

W(r)  = 

w^(r-l) 

W2(r-1) 

..  Wp(r-1) 

Wj(r-p+2) 

W2(r-p+2) 

. .  w  (r-p+2) 

P  J 

V(r),  y(r)  being  the  same  in  form  as  W(r)  with  replaced  by 
respectively. 

Now  each  of  the  sets  of  functions  w^(r),  v^(r),  y^(r)  are  linearly 
independent,  so  that  there  can  be  at  most  a  finite  number  of  values  of 
X  for  which  any  of  the  matrices  W(r),  V(r)  or  Y(r)  are  singular.  There 
is  no  real  loss  of  generality  if  we  suppose  that  W(r),  V(r),  Y(r)  are 
all  non-singular  for  all  r  >  p  -  1,  for  otherwise  we  simply  ignore  a 
finite  number  of  terms  at  the  beginning  of  the  series,  convergence  not 
being  affected.  Similarly,  there  is  no  loss  of  generality  in  assuming 
that  (2.2)  holds  for  all  r  >  p  -  1. 

This  means  that  not  only  do  we  have  (as  we  know  already)  the 
constant  vector  A  such  that 

U(r)  =  W(r)  A  ,  (3. 6) 

but  also  uniquely  determined  vectors  B(r),  C(r)  given  by 

U(r)  =  V(r)B(r)  =  Y(r)C(r)  .  (3.7) 


Now  we  have,  clearly,  as  r  «  , 
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W(r)  ~  V(r)  , 

so  that 

B(r)  =  [V(r)l"^U(r)  =  [V(r)]"^W(r)A  ~  A  .  (3.8) 

Moreover,  if  we  Introduce  the  diagonal  matrix  D(r)  given  by 

<^1  <r  <r 

D(r)  =  <r  ',r  ^  ,r  P  >  ,  (3.9) 

'"i  ‘"i 

then  from  the  fact  that  as  r  -♦  <»  ,  for  fixed  k  ,  (r  +  k)  r  we 
have  immediately 

V(r)  -  D(r)Y(r)  , 

and  hence 

C(r)  =  [Y(r)]"^U(r)  ~  D(r)[V(r)]"^U(r)  , 

=  D(r)[V(r)]“^W(r)  A  , 

~D(r)A  .  (3.11) 

Now  by  means  of  (3.2),  (3.4)  we  reformulate  equation  (1.2)  as 

U(r)  =  K(r)U(r  -  1)  ,  (r  >  p  -  1)  .  (3. 12) 

the  "starting  values"  u(0),  ...,  u(p-l)  constituting  the  "starting  vector" 

U  (p  -  2) .  Hence  we  have 

U(r  +  1)  =  K(r  +  l)K(r)U(r  -  1)  , 
etc. ,  and  if  we  write,  for  brevity, 

L(n,s)  =  K(n+ s)K(n+ s  -  l)...K(n+  1)  ,  (s  >  1),  L(n,0)  =  1  ,(3.13) 
we  have 


U(n  +  s)  =  L(p  -  2,r  -  p  +  2)U(p  -  2) 


(3. 15) 
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which,  by  expressing  U(r)  in  terms  of  the  known  elements  of  the 
matrices  K(r)  and  the  elements  of  the  starting  vector  U(p  -  2),  provides 
a  formal  solution  of  the  equation  (3. 12).  It  tells  us  nothing,  however, 
about  the  convergence  of  the  corresponding  series  S  uWx*^  .  So  we 
introduce  the  vector  B(r)  from  (3.  7)  and  have 

B(r)  =  [V(r)]"^L(p  -  2,r  -  p  +  2)U(p  -  2)  ,  (3.  16) 

and  (3.  8)  then  gives 

A  =  lim  [V(r)]"^L(p  -  2,r  -  p  +  2)U(p  -  2)  .  (3. 17) 

r-*  00 

So  the  condition  a^  =  0  is  equivalent  to  the  condition  that  the  l.^*'  element 
in  [V(r)]”^L(p  -  2, r  -  p  +  2)U(p  -  2)  should  tend  to  zero  as  r  oo.  To 
express  this  in  matrix  form,  let  denote  a  row  vector  of  £  elements, 
in  which  the  n  element  is  unity  and  the  rest  all  zero,  e.g. 

Jj  =  [l,  0,  0, . . . ,  O]  ,  ^2  ~  •  •  • »  (3. 18) 

then  we  have 

a^  =  A  ,  (3. 19) 

so  the  condition  a^  =  0  is  equivalent  to 

lim  J  [V(r)]"^L(p  -  2,  r  -  p  +  2)U(p  -  2)  =  0  (3. 20) 

r  -*  00 

In  this  formula,  it  should  be  noted,  we  are  net  using  the  full  solutions 
w^(r)  but  only  their  dominant  terms  v^(r)  which  are  easily  found. 

The  radius  of  convergence  of  2  u(r)x^  will  thus  be  in  general. 
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but  will  be  if  the  equation 

lim  Jj[V(r) ]  "^L(p  -  2,  r  -  p  +  2)U(p  -  2)  =  0  ,  (3. 21a) 

r  -♦00 

is  satisfied,  and  if  the  condition 

lim  J2[V(r)]  '^Kp  -  2,r  -  p  +  2)U(p  -  2)  =  0  ,  (3.  21b) 

r  -*■«) 

is  also  satisfied,  and  so  on. 

The  simplest  non-trivial  example  of  this  is  that  in  which  p  =  2, 

Pi  <  P2  ,  and  the  coefficients  a^fr),  a^{r)  contain  a  single  parameter  \, 
say,  the  ratio  being  given  either  explicitly  or  in  terms  of  X. 

(This  is  the  situation  which  arises  in  the  study  of  Mathieu  functions, 
spheroidal  wave  functions,  and  Lame  functions,  hitherto  studied  by  the 
continued-fraction  technique).  The  problem  is  then  to  determine  X  so 
that  the  radius  of  convergence  may  be  augmented  from  Pj  to  .  In 
such  circumstances  (3.21a)  provides  an  equation  for  calculating  the 
possible  values  of  X  and,  in  the  particular  cases  mentioned,  is 
generally  known  as  the  characteristic  equation.  To  keep  in  line  with 
this  practice,  we  shall  call  (3. 20)  the  i^^  characteristic  equation  ; 
for  brevity,  denote  it  by  =  0.  Then,  reverting  to  general  values  of  £  , 
suppose  there  are  (p  -  1)  parameters  X  X  . . . ,  X  involved 

in  the  problem.  The  radius  of  convergence  will  be  p^  in  general;  in 
order  that  it  should  be  augmented  to  p^  the  one  relation  =  0  must 
hold  among  the  parameters  X  ;  if  the  further  relation  E2  ”  ® 
holds  the  radius  of  convergence  is  augmented  to  p^,  etc. ;  ultimately 
the  radius  of  convergence  can  be  increased  to 


if  p  -  1  relations 
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hold  between  the  p  -  1  parameters. 

A  further  form  of  the  characteristic  equations  can  be  obtained,  in 
which  we  do  not  even  have  to  form  the  matrix  V(r)  but  need  only  the 
still  simpler  matrix  Y(r).  By  similar  reasoning  to  that  which  led  to 
(3. 16)  we  have 


C(r)  =  [Y(r)]"^L(p  -  2,r  -  p  +  2)U(p  -  2)  ,  (3.  22) 


from  which  the  coefficients  c^W  can  be  obtained  immediately  (in 

matrix  notation,  we  have  c^(r)  =  J^C(r)).  Now  we  know  that  as 

""i 

r  -►  00,  c^(r)~  r  a^  ,  so  the  condition  c^(r)-»  0  does  not  necessarily 

imply  a^  =  0.  However, 

,  ,v,l/r  ^/r  1/r 

{c^(r)  }  ~  r  a^  ' 

1  if  a.  #  0  , 


0  if  a^  =  0  , 


so  that 


th 


>l/r 


{c^Wr'*  -0  a^  =  0  , 


and  the^  characteristic  equation  can  be  written 


(3. 23) 


lim  {Jj[Y(r)]"^L(p  -  2,r  -  p  +  2)U(p  -  2) }  =  0 

r-^00 


(3.  24) 


Finally,  in  the  specially  favourable  case  when  all  the  are 
different,  yet  another  form  of  the  characteristic  equation  is  available. 
For  in  such  a  case  we  have,  if  a^  #  0  , 
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u(r)  ajWj(r)  =  Cj(r)yj(r)  ; 


(3.  25) 


hence 


u(r+Al 

u(r) 


Cj(r+1)  yj(r+l) 

Cj(r)  yj(r) 


M»l 

but  yj(r  +  l)/yj(r)  ^  r  t^e 
-1*1'^,  -1  -“l  u(r+l) 

*1  «  “W) 


,  SO 

aj(r+l) 


(1  +  ^  ~  1  .  (3. 26) 


If,  however,  a^  =  0,  a^  ^  0,  u(r)  ^  and  we  find  by  similar 

working  that 


1  u(r+l) 
u(r) 


(t2/tj)e 


9 


and  the  right  hand  side  of  this  approaches  the  limit  0  if  M’2  > 
^*2^^1^  if  |ij  =  (12*  Hence  the  condition  a^  =  0,  a2  #  0 


lim 
r-^  00 


**1*^  -1  **1  u(r+l)  ^  r 

r  tj  e  , ,  =  0  if  |Xj>  ^3  , 


(3.  27) 


Since  u(r  +  1)  =  JjU(r),  u(r)  =  J2U(r),  and  U(r)  is  given  by  (3.15),  this 
provides  an  alternative  form  of  the  first  characteristic  equation. 


4.  Direct  determination  of  the  possible  radii  of  convergence 


From  the  results  in  §  3  we  obtain  easily  the  following  working  rule 
for  determining  the  possible  radii  of  convergence. 

In  the  relation  (1.2),  divide  by  the  highest  power  of  jc  occurring, 
then  let  r  -*oo,  and  let  the  limiting  form  of  (1.2)  be 

wu(r  +  1)  =  a,u(r)  +  a,u(r  -  1)  +  . . .  +  c  u(r  -  p  +  1)  ,  (4.1) 

12  p  ' 


u,  Oj,  Q^, . . .,  ttp  being  constants,  not  all  zero. 

Then  the  possible  radii  of  convergence,  p^,  are  the  roots  of  the 
equation 

-  u  +  a,  p  +  +  . . .  +  O'  p^  ^  =  0  ;  (4.  2) 


this  equation  is  to  be  regarded  as  essentially  of  degree  (p-1),  so  that 
infinite  roots  (if  they  occur)  must  be  counted. 

The  justification  of  this  rule  depends  on  the  fact  that  the  existence 
of  a  solution  of  (1.  2)  of  the  form  (1.4)  or  (1.  5)  with  pi  >  0  implies  the 
existence  of  a  zero  root  of  (4.  2),  the  existence  of  a  solution  with  p.  <  0 
implies  an  infinite  root  of  (4.  2),  while  for  solutions  with  p  =  0  the 
possible  values  of  ^  are  given  by  the  reciprocals  of  the  finite  non-zero 
roots  of  (4. 2),  as  one  easily  sees  by  actual  substitution  of  the  trial 
solutions  (1.4)  or  (1.5). 

Hence  the  possible  radii  of  convergence  of  the  series  Z  u(r)x'^ 
can  be  written  down  at  once  from  the  recurrence  relation  (1. 2);  this 
generalizes  the  rule  given  by  Perron  in  (9)  Vol.  II,  §  46. 


f 
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5.  Determination  of  the  coefficients 

One  might  think  that  when  the  starting-values  or  the  parameters 
have  been  determined  in  such  a  way  as  to  secure  the  desired  radius  of 
convergence,  the  problem  would  be  solved  in  the  sense  that  the 
coefficients  u(r)  can  be  found  successively  by  use  of  (1.2).  However, 
in  the  simplest  non-trivial  case,  namely  p  =  2,  it  has  long  been  known 
that  this  is  not  so  (10)>  Unless  we  achieve  the  impossible  by  determining 
the  exact  values  of  any  parameters  involved,  and  making  no  approxi¬ 
mations  whatever  at  any  stage,  the  small  errors  involved  will  Inevitably 
buildup,  re-introducing  undesired  functions  Wj(r),  W2(r),  etc.  into 
u(r) ,  and  so  destroying  the  augmented  convergence. 

To  avoid  this,  we  need  to  have  expressions  which  give  the  U(r) 
(and  hence  the  u(r))  directly.  The  process  of  obtaining  these  leads 
to  yet  another  form  for  the  characteristic  equations. 

For  simplicity,  suppose  first  that  a  (r)  #  0  for  all  r  ,  so  that 

P  ” 

the  matrix  K(r)  is  non-singular  for  all  x  I  the  exceptional  case  will 
be  dealt  with  in  the  next  section. 

Suppose  we  wish  our  u(r)  to  be  a  predetermined  combination  of  the 
Wj(r)  -  that  is  to  say,  the  a^  are  to  have  pre-assigned  values.  (For 
instance,  we  may  wish  to  have  a^  =  0  for  i  =  1, 2, . . . ,  p  -  1  in  order  to 
secure  maximal  convergence  for  Z  u(r)x*^. )  Now  we  have,  from  (3.7), 

(3.  14) 


(5.  1) 


L(n,  s)U(n)  =  V(n  +  s)B(n  +  s) 


f 
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hence 

B(n  +  s)  =  [V(n  +  s)]"^L(n,  s)U(n)  .  (5.2) 

Now  let  s  -*  oo;  then  by  (3.8)  the  left  hand  side  tends  to  the  definite 

limit  A,  so  that  the  limit  of  the  right-hand  side  certainly  exists.  At 

present  we  are  excluding  the  case  a  (r)  =  0,  so  all  the  K(r)  are 

P 

non-singular,  L(n,  s)  is  also  non- singular  and  hence 

lim  [L(n,  s)]  ^V(n  +  s) 

s  -►  OO 

exists.  But  from  (3.4)  we  have 

U(n)  =  [L(n,  s)]  ^V(n  +  s)B(n  +  s)  (5.  3) 

so  that  on  letting  s  -►  oo  we  have  (since  B(n  +  s)  -►  A), 

U(n)  =  lim  [L(n,  s)]  V(n  +  s)A  .  (5.4) 

S  -*  00 

In  passing,  it  may  be  noted  that  we  can  easily  verify  directly  the 
fact  that  U(n),  as  given  by  (5.4),  does  satisfy  the  difference  equation 
in  its  vector  form  (3.12).  Observing  that 

L(n,  s)  =  L(n  +  1,  s  -  l)K(n  +  l)  , 

which  is  an  immediate  consequence  of  the  definition  of  L(n,  s)  in  (3. 13) , 
we  obtain 

U(n)  =[K(n  +  1)]’^  Urn  [L(n  +  1,  s  -  l)]"^V(n  +  s)A 

s  -►  OO 

=  [K(n  +  l)]'V(n  +  1)  , 

so  that  K(n  +  l)U(n)  =  U(n  +  1),’  which  is  the  same  as  (3. 12). 
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Going  back  to  relation  (5.4),  we  bring  in  another  expression  for 
U(n),  namely  that  obtained  in  (3.  15),  and  combining  the  two  we  have 

L(p  -  2,  n  -  p  +  2)U(p  -  2)  =  lim  [L(n,  s)]"  V(n  +  s)A  .  (5.  5) 

S-*  00 

In  this,  n  may  be  any  integer  >  p  -  1  ;  for  n  =  p  -  1  we  have  a  simpler 
form 

U(p  -  2)  =  11m  [L(p  -  2,  s)]  V(p  -  2  +  s)A  .  (5.6) 

Let  us  consider  the  significance  of  this  equation.  Each  side  is 
a  £  -  element  column  vector,  so  that  (5.  5)  or  (5.  6)  contains  effectively 
£  linear  equations.  It  provides  the  necessary  conditions  which  must 
be  satisfied  in  order  that  the  elements  a^  of  the  vector  A  shall  have 
pre-determined  values,  these  conditions  involving  the  £  elements  of 
the  starting  vector  U(p  -  2)  and  the  elemencs  of  the  matrices  K(r)  which, 
as  remarked  already,  may  contain  parameters.  If  we  wish  convergence 
to  be  augmented  from  to  then  we  must  have  a^  =  0,  but  the 
remaining  a^  are  arbitrary,  so  that  (5.  5)  will  provide  precisely  one 
condition  to  be  satisfied,  which  is  thus  equivalent  to  the  first  character¬ 
istic  equation  (3.  21a).  If  we  wish  convergence  to  be  augmented  to  p^, 
we  must  have  a^  =  a^  =  0,  so  that  (5.  5)  then  has  in  it  two  conditions, 
equivalent  to  the  first  and  second  characteristic  equations  and  so  on. 

For  maximal  convergence,  we  must  have  a^  =  0  for  1  <  p,  so  that 
(5.  5)  provides  (p  -  1)  conditions. 


6.  Case  of  a  singular  K(r) 


We  have  now  to  show  how  the  analysis  of  §  5  is  affected  when 

QfpCr)  vanishes  for  some  x  >  causing  the  corresponding  matrix  K(r) 

to  become  singular.  For  simplicity,  we  shall  deal  only  with  the  case 

when  there  is  a  single  value  of  r  ,  say  r  =  N,  for  which  a  (r)  =  0. 

“  P 

It  is  still  possible,  then,  to  construct  a  solution  having  prescribed 
values  of  the  a^,  by  using  (5.  4)  with  n  >  N.  For  n  <  N  we  must  use 
the  equivalent  of  (3.  15),  namely 

U(n)  =  L(p  -  2,  n  -  p  +  2)U(p  -  2)  (6. 1) 

Since  this  is  used  only  for  the  calculation  of  a  finite  number  of  U(n), 
the  small  errors  inevitably  introduced  will  not  build  up  to  destroy 
augmented  convergence.  The  equation  (5.  5)  still  holds  good  provided 
we  have  n  >  N,  and  can  be  used  in  the  same  way  to  obtain  the  various 
characteristic  equations. 

One  other  possibility,  of  considerable  interest,  arises  in  this  case 
but  not  in  that  of  §  5.  We  have,  for  the  column  vector  U(  N  -  1)  , 

U(N  -  1)  s{u(N),  u(N  -  1), . . . . ,  u(N  -  p  +  1)} 

with 

U(N  -  1)  =  L(p  -  2,  N  -  p  +  l)U(p  -  2)  (6.  2) 

Suppose  we  choose  the  (p  -  1)  parameters  in  the  problem  so  as  to  make 

u(N)  =  u(N  -  1)  =  . . . .  =  u(N  -  p  +  2)  =  0:  then  since  a  (N)  =  0  we 

P 

have  u(N  +  1)  =  0  also,  and  hence  U(N)  =0.  But  we  have 
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U(N  +  1)  =  K(N  +  l)U(N),  etc.,  and  consequently  U(n)  =  0  for  all 
n  >  N.  The  series  (1. 1)  then  reduces  to  a  polynomial  of  degree  (N  -  p  +  1) 
in  X.  Since  a  polynomial  is,  of  course,  convergent  everywhere  this  may, 
in  certain  circumstances,  provide  a  means  of  augmenting  the  convergence 
from  Pp  (if  this  is  finite)  to  infinity. 

As  an  illustration  of  this,  consider  the  spheroidal  wave  equation 

(1  -  x^)  -  2x  +  {  X.  -  ti^(l  -  x^)  ^  +  Y^(l  -  x^)}  w  =  0  ;  (6.  3) 

dx 

2 

regarding  y  as  fixed  and  \,  pi  as  disposable  parameters.  Assuming 
a  solution  of  the  form 

oo 

S  u(r)x'^  (6. 4) 

0 

the  general  recurrence  relation  is  found  to  contain  foxir  terms  (1.  e.  p  =  3) 
and  to  be 

(2r+  l)(2r+2)u(r+  1)  =  (8r^  +  pi^  -  \  -  Y^)u(r)  + 

+  {\  +  2y^  -  (2r  -  l)(2r  -  2)}  u(r  -  1)  +  Y^u(r-  2)  . 

2 

Dividing  through  by  4r  and  letting  r  -»  «>  we  obtain,  in  the  notation  of 
S4, 

w  =  1,  =  2,  =  -  1,  =  0  ,  (6.6) 

giving  Pj  =  P2  “  P3  =  “•  Hence  by  suitable  choice  of  \  and  pi  we 

can  obtain  an  infinite  series  solution  which  is  an  integral  function  of  x  . 

2 

Here  a  (r)  =  y  ^0  so  the  general  method  of  §  5  applies. 

P 

Now,  however,  consider  the  degenerate  case  where  y  =  0  (!•  the 
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Associated  Legendre  equation).  The  situation  is  now  different,  since 

we  have  p  =  2,  and  the  possible  radii  of  convergence  are  1,  1  so  that 

a  solution  which  is  an  infinite  series  of  the  form  (6.4)  cannot  possibly 

converge  outside  |x|  =1.  The  only  way  to  obtain  a  solution  with 

better  convergence  properties  is  to  make  use  of  the  fact  that  now, 

since  p  =  2,  we  have  a  (r)  =  \  -  (2r  -  l)(2r  -  2),  so  that  if  we  choose 

P 

X.  as  (2N  -  1)(2N  -  2),  «  (N)  =  0  and  we  have  the  situation  of  this 

P 

paragraph.  By  suitable  choice  of  |x  we  can  make  u(N)  vanish,  and 

2 

the  solution  then  becomes  a  polynomial  of  degree  (N  -  1)  in  x  . 
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